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Abstract 

The kinetics of phase and population relaxation in quantum systems induced by noise with anoma- 
lously slowly decaying correlation function P{t) ~ (wt)'" , where < a < 1 is analyzed within con- 
tinuous time random walk approach (CTRWA). The relaxation kinetics is shown to be anomalously 
slow. Moreover for q < 1 in the limit of short characteristic time of fluctuations w"^ the kinetics 
is independent of w. As q ^ 1 the relaxation regime changes from the static limit to fluctuation 
narrowing. Simple analytical expressions are obtained describing the specific features of the kinetics. 

Introduction. The noise induced relaxation in quantum systems is very important process observed 
in: magnetic resonance PP, quantum optics and nonlinear spectroscopy 0, etc. 

Often these processes are analyzed assuming conventional stochastic properties of the noise: fast 
decay of correlation functions and short correlation time Tc . In the absence of memory the relaxation 
is described by very popular Bloch-type equations. As for memory effects, they are also discussed (within 
the Zwanzig projection operator approach jS]), however, in the lowest orders in the fluctuating interaction 
V inducing relaxation 

Recently strong attention has been drawn to the processes governed by noises with anomalously slowly 
decaying correlation functions P{t) ^ t^" with a < 1. They are discussed in relation to spectroscopic 
studies of quantum dots references therein). Similar problems are analyzed in the theory of 

stochastic resonances j^. 

Such anomalous processes cannot be properly described by methods based on expansion in powers of 
V. The goal of this work is to analyze the corresponding anomalous relaxation within the continuous 
time random walk approach (CTRWA) |H| with the use of the recently derived non-Markovian stochastic 
Liouville equation (SLE) |^ which enables one to describe relaxation kinetics without expansion in V. In 
some physically reasonable models it allows for describing the phase and population relaxation kinetics 
in analytical form even for multilevel systems. In particular, the kinetics is shown to be strongly non- 
exponential . 

General formulation. We consider noise induced relaxation in the quantum system whose dynam- 
ical evolution is governed by the hamiltonian 

H{t)=H, + V{t), (1) 
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where Hg is the term independent of time and V{t) is the fluctuating interaction, which models effects of 
the noise. The evolution is described by the density matrix p{t) satisfying the Liouville equation [h = 1) 

p ^ -iH{t)p, with Hp = [H, p] = [Hp - pH]. (2) 

y(t)-fluctuations are assumed to be symmetric {{V) = 0) and result from stochastic jumps between 
the states \x,^) in the (discrete or contimmm) space {x^} = {x} with different V = V^, and H — H^ (i.e. 
different V = V^ = [K, . . .] and H = H^): 

V = y \x^)%{xu\ and7T: = V \x^)H^{x^\. (3) 

Hereafter we will apply the bra-ket notation: |fc) and \kk') = \k){k'\ for the eigenstates of H (in the 
original space) and H (in the Liouville space), respectively, as well as notation \x) for states in {x}-space. 

The macroscopic evolution of the system under study is determined by the evolution operator TZ{t) 
in the Liouville space averaged over T^(t)-fluctuations: 

p{t) = n{t)p, withn{t) = V g{x,x,\t)Peix,), (4) 

where Q{x,x'\t) is the average evolution operator and Pe{x) is the equilibrium distribution in {a;}-space. 

Non-Markovian l/(i)-fluctuations will be described by the CTRWA [which leads to the non-Markovian 
SLE in] for g{t)]. It treats fluctuations as a sequence of sudden changes of V. The onset of any particular 
change of number j is described by the matrix Pj-i (in {a;}-space) of probabilities not to have any 
change during time t and its derivative Wj-i(t) = —dPj-i{t) / dt. These matrices are diagonal and 
independent of Pj-i{t) ^ P(t), Wj_i(t) = W(t) = ~dP{t)/dt, {j > 1) , except Fo(0 = Pz{t) and 
Wo{t) = Wi{t) — —dPi{t)/dt depending on the problem considered. For non-stationary (n) and stationary 
(s) fluctuations 8 W^{t) = Wjt) = W{t) and W^{t) = Ws{t) = f-^ J^°°dTW{T), respectively, where 
iw = f^d,TTW{T) is the matrix of average times of waiting for the change [5]. 

In what follows we will mainly operate with the Laplace transforms denoted as Z{e) — f^^dt Z{t)e~^* 
for any function Z{t). In particular, noteworthy is the relation Pj{(-) = [1 — Wj{e)\/t and suitable 
representations 

W{e) = [1 + $(e)]-i and P(e) = [e + e/^e)]~^ (5) 

in terms of a diagonal matrix $(e) with $(e) ~ (e/i&)", where u) is a constant matrix and a <1 (see 
below). 

Evolution in {a;}-space is governed by the jump operator C = \ — V\n which V is the the non-diagonal 
matrix of jump probabilities. Evolution leads to the equilibrium state le^;), satisfying eq. C'w°'\ex) = 
and represented as \ex) = '}2xPe{x)\x) with (e^r] = '}2xi^\ 0- Note that [see eq. Q] 

n{t) - {ex\Q\ex) = {G). (6) 
The CTRWA leads to the non-Markovian SLE for g{x,x^\t) 9 . Solution of this SLE yields 

g = p,in) + f2-i<l(i7)[$(f2) + c]-^vw,in), (7) 

where 

C = 1 -V a.nd n = e + iH. (8) 
In particular, in the case of n-fluctuations {Wi — W) 

g = gn^n-'H^ + c)-\ (9) 
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For s-fluctuations {W^ =Ws)G = (l^^ - QntiClL)-^. 

Hereafter, for brevity, we omit the argument VL of all Laplace transforms if this docs not result in 
confusions. 

Useful models and approaches. 

Sudden relaxation model (SRM). The SRM ^ assumes sudden equilibration in {xj-space described 
by operator 

^=(l-|eo)(eo|)g'\ Q^l~Y..Po{x)\x){xl (10) 
in which |eo) = '^xPo{x)\x) , (eo| — '^x{x\- For this C 



le^;) = gleo) with q = " / {eo\Qw "|eo) 
and {ex\ = (eo|. In the model (fTn|l one gets for any Wi 



(11) 



(12) 



where Pq^ = (1 - Wq,)ICi, Wq^{1 + $Q)-i, and Wq^ = Wi{WQ/W). 

Short correlation time limit (SCTL). In practical applications of special importance is the SCTL for 
F(i)-fluctuations in which eq. (|12|) can be markedly simplified. It corresponds to large characteristic 
rates Wc of the dependence $(17) = ^(Cl/wc)'- Wc ^ ll^ll- In this limit the relaxation kinetics is described 
by the first terms of the expansion of $(i7/wc) in small Cl/wc, since <I*(e) is the increasing function of e 
with $(e) "^-^ . Some important general conclusions, however, can be made independently of the form 
of $(12) (see below). 

Models for quantum evolution and fluctuations. The obtained general results are conveniently illus- 
trated with the quantum two-level model and the stochastic two-state SRM for F(t)-fluctuations. 

Quantum evolution of the two-level system is governed by the hamiltonian (assumed to be real matrix) 
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Vn -Vd 



(13) 



The two-state SRM suggests that fluctuations result from jumps between two states (in {a;}-space), 
say, |a;+) and whose kinetics is described by 



£ = 2(1 - \ex){ex\) with \ex) = ^||a;+) + \x-)). 
Below we will consider two examples of these models: 

1) Diagonal noise [IT)] : ujg = 0, Vn — 0, Vd — loo{\x+){x+ \ — \x-){x-\), and 

H.=± = ±^M\+){+\ l-X-l); 

2) Non-diagonal noise: Vd = and V„ = v{\x^){x+\ — |a;_)(a;_|), so that 

H,^±=Hs ±v{\+){-\ + |-)(+|). 



(14) 



(15) 



(16) 



The first model describes dephasing while the second is useful for studying population relaxation. 
In the model H13I) dephasing and population relaxation are characterized by two functions: the spec- 
trum I{uj) and the difference of level populations N{t): 



^ iRe(s|7^(iw)|s) and N{t) = (r^|7^(^)|r^) 



(17) 
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with |,s) = ^11 + -) + 1-+)) and \n) = ^|| + +> - |--)). 

General results in the SCTL. Within the SCTL (||T^||/?i;c 1) especially simple results are 
obtained for ||_ffs||/?i;c ^ 1- In the lowest order in ||$(ii/wc)|| ^ 1 

n^Un^ {q-^Q^^^{n))/{q-^Q^{n)). (is) 

This formula holds for any initial matrix Wi and, in particular, for s- fluctuations, if \\tw\\ ~ l/u^c 
l/||f2||. 

The more complicated SCTL-case ||iJs|j/ii;c — 1 can be analyzed by expanding Q in powers of the 
parameter ^ = ||F||/||iJs|| ^ 1- In particular, within the general two-level model [eq. (|13ll ] with Vd = in 
the second order in ^ the diagonal and non-diagonal elements of p{t) are decoupled and the corresponding 
elements of TZit) are expressed in terms of the universal function 

1 1"^°° e**"* 



Rk(t) = — / de : (19) 

{fim)\fi) = e-^^-'Rk,{t), (/. = n,+-, -+), (20) 
where — kn — 2Re(fc^ ), and 

k+- = fc*+ = ^iO-^Vnq-'[l - WQ{2lLU,)]Vn). (21) 

Anomalous fluctuations. The simplest model for anomalous fluctuations can be written as JI] 

$(e) = (e/wr, (0 < a < 1), (22) 

where w is the matrix of fluctuation rates diagonal in |a:;)-basis. For simplicity, w is assumed to be 
independent of x, i.e. w = w(= Wc)- The model H22|l describes anomalously slow decay of the matrix 
W{t) ^ l/f^+" (very long memory effects in the system JT,), for which only the case of n-fluctuations is 
physically sensible. 

In the SCTL (fH^ the model yields the expression 

7^„(e) = (A"-i(e))(Sl"(e))-i with n{e) ^e + iH (23) 

which shows that TZn{^) [and 7^„(i)] is independent of the characteristic rate w. For a = and a = 1 eq. 
H23|l reproduces the static and fluctuation narrowing limits 1 : 7^,i(e) = (rj^^(e)) and TZn{e) = l/(r2(e)), 
respectively. 

Of certain interest is the limit a — * 1 in which formula H23|l predicts the Bloch-type exponential 
relaxation: 

7^„(e) « [e + iH, + ia~l){nHn)),^o]-\ (24) 

controlled by the relaxation rate matrix Wr — {a — l)Re(f2 ln(17))e^0: and accompanied by frequency 
shifts represented by h ~ i{a — l)Im(f2 ln(r2))e_>o- However, the matrices Wr and h (unlike those in the 
conventional Bloch equation) arc independent of the characteristic rate w of l/(i)-fluctuations. 

Dephasing for diagonal noise. In the model ()15|> the spectrum I{uj) can be obtained in the general 
SRM (dnj: 

""(^a)2 + (^a)2 + 2V'2;V'+COs(7ra)' ^ ' 
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where V'±('^) = {\^~'^Vd\^d[±{uj—2Vd)]) with ^(z) being the Heaviside step-function and = sin(7ra)/7r. 

In the two-state SRM ^ I{lo) = \naL0Q'^9{y) {y + y^^ + 2)/[y" + y^" + 2cos(7ra)], where y = 
{uq + u)) / {ljJq — Lj) (see also ref. According to this formula anomalous dephasing (unUke conventional 
one pP) leads to broadening of I{lo) only in the region < ojq and singular behavior of I{uj) at — > ±wo: 
l/{uj ± wq)"'^ For a > Qfc ~ 0.59 [ac satisfies the relation ac = cos(7rac/2)] the two-state- 
SRM formula also predicts the occurrence of the central peak (at uj ~ 0) 'W Lorenzian shape and 
width WL ~ Wo cos(7ra/2)/-\/a2 — cos2(7rQ:/2): I{uj) « tan(7ra/2)(jjQ"^/[l + (uj/wl)'^], whose intensity 
increases with the increase of a — ac (Fig. 1). At a ~ 1 the parameters of this peak are reproduced 
by eq. H24|l in which (r21n(r2)}e = —{tt/2)u!o. The origination of the peak indicates the transition from 
the static broadening at a <C 1 to the narrowing one at a ^ 1 [see eq. (jSSJ]- For systems with complex 
spectra this transition can, of course, be strongly smoothed. 

Dephasing for non-diagonal noise. The model allows one to reveal some additional specific 
features of dephasing. We restrict ourselves to the analysis of the case \\Hs\\ ^ uig ^ w and the 
most interesting part of the spectrum at ^ ujs- Equations H19(l and 1)20(1 show that in this case 
the elements (/i|7^(t)|/i), (fi = H — , — h), which describe phase relaxation are given by {fi\TZ{t)\fi) — 
e-^'^^^tEa[-k^{wt)°'], where Ea{-z) = {2TTi)-^ f_^^dy / {y + zy^-°') is the Mittag-Lefflcr function 
Therefore for |a;| ^ tOs 

I{uj) ^ Ia{ujs + i^) + Iq{^^s ~ (^), (26) 

where 

/o(w) = no sm<j)^{\x\^+°' + + 2\x\ coscj)^)-^ (27) 

with X — w/(|fc_| \'^/°'w), no = (7r|fc_| \^^°''w)~^, and 



7ra . , , 
bx = — — h sign(a:) arctan 



sin(i7ra) 



(28) 



_cos(i7rQ!) + 2 " LOs/w_ 

Formula ((26(1 predicts singular behavior of I{uj) at a; ^ ±ojs- ^(w) ^ l/\uj ±ujs\^~" , and slow decrease of 
I{uj) with the increase of | w ± Ws | : I{uj) l/\uj ztojsl^^"- 

In the limit Us/w <^ 1 (j)x ~ ■Ka6{x) so that Io{i-i->) ~ ^(^)- This means that for ajs/w <^ 1 the 
spectrum /(w) is localized in the region |a;| < lUs and looks similar to I{lu) for diagonal dephasing at 
a < ac (see Fig. 1). For aJs/w > 1, however, I{lu) is non-zero outside this region, moreover, in the limit 
uJs/w ^ 1 the spectrum Iq{uj) becomes symmetric: loi^^) = loi^^), similar to the conventional spectra. 

It is also worth noting that for ujs/w <C 1 functions {fi\TZ{t)\fi) and I{w) are independent of w [in 
agreement with eq. ((18(l ] since fc^ ~ (us/w)" and kf^iwt)" ^ (w^t)". In the opposite limit, however, 
kn ~ so that the characteristic relaxation time ^ w^^. 



Population relzixation. Specific features of anomalous population relaxation can be analyzed with 
the model of non-diagonal noise 1(16(1 . 

In particular, in the limits ^ ujg >w and 1 — a ^ 1 with the use of eqs. ((T^ . lf^ and ((21jl one 

gets 

N{t) = Ea[-kn{wt)"'] and N{t) = e-'"^*, (29) 

respectively, where Ea{—x) is the Mittag-Leffer function defined above and Wa ~ fc„(Q ^ l)w 1 — a. 
The first of these formulas predicts very slow population relaxation at t > Tr = w~^(fc„/ui)^/": N{t) ~ 
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Similar to I{lu) the function N(t) is, in fact, independent of w in the hmit lUs/w <^ 1 because in 
this case kn ^ {lOs/w)"'. In the opposite hmit los/w > 1 the characteristic time population relaxation is 
~ since k„ is independent of w as with the phase relaxation. 
In the limit ||i?s||7 «C w one obtains 



where np{e) = [{e + 2iEo)'^ + {e-2iEof]/2 and Eq = ^v^ + lo^ /A. Naturally, in the corresponding limits 
expression (|30l) reproduces formulas H29|) with fc„ « 2"^^ cos(7ra/2)(_Eo/w)" and Wa ~ 7r(l — a)u^/i?o 
(Fig. 2). Outside these limits N{t) can be evaluated numerically (some results are shown in Fig. 2). In 
general, N{t) is the oscillating function (of frequency ~ Eq) with slowly decreasing average value and 
oscillation amplitude: for E^t ^ 1 N{t) ^ (except the limit a 1). 

Concluding remarks. The presented analysis of relaxation kinetics in quantum systems induced by 
anomalous noise demonstrates a number of peculiarities of this kinetics. The peculiarities are analyzed 
with the use of the two-level quantum model, as an example, though the observed anomalous effects can 
show themselves in more complicated multi-level quantum systems. The proposed theoretical method is 
quite suitable for the analysis of these systems. This work is currently in progress. 

Noteworthy is that in some limits the theory developed predicts relaxation kinetics described by the 
Mittag-Lefher function (wt)"]. Following a number of recent works (for review see ref. 11 ) this 
kinetics can be thought as a result of the anomalous Bloch equation with fractional derivative in time. 
For brevity we have not discussed the corresponding representations. 

It is also interesting to note that with the increase of a the effects of anomaly of fluctuations decrease 
but still persist. To clarify them let us briefly consider the model $(e) — (e/w) + C(e/w)^+", in which 
< a < 1, and w and C are the constants with C 1 [small value of ( ensures that W{t) > 0]. Possible 
effects can be analyzed within the SCTL with the use of eqs. (|18|) - (I20|) . For example, in the limit 
< 1 one obtains formula U w [e + iH^ + (w^"' {{in)^+°' - {iHs)'^+°')]~'^ predicting the Bloch- 
type relaxation of both phase and population, but with the rate Wr — Cw^"Re((i7i)^+" — [iHsY^") 
which depends on w as i.e. slower than in the conventional Bloch equation [Wr ~ 1/w 1 ). 

Analysis also shows that in the expression for TZ the terms ^ w(e/u')^+" occur as well. They lead to 
the inverse power-type asymptotic behavior of (/i|i?(i)|/i) ~ l/t^+" observed, however, only at very long 



In our brief analysis we neglected the effect of possible natural width of lines corresponding to the 
additional slow exponential relaxation in the system. It is clear that the developed method allows one 
take in account these effects straightforwardly in case of need. 
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(30) 



times t^w 



-1 
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Figures. 

Fig. l.The spectrum I{x) = I{lo)ujq, where x = lo/ujo, calculated in the model H15|) [using eq. (|25|l ] 
for different values of a: (1) a = 0.5, (2) a = 0.7, (3) a = 0.8, and (4) a = 0.9. 

Fig. 2. Population relaxation kinetics N{t), where r = Eot, calculated with eq. H30|) (a) for large 
a and different r = 2t;/cj,: (l)-a = 0.95, r = 1.0; (2)-a = 0.95, r = 2.0; (3)-a = 0.88, r = 1.0; 
a — 0.88, r — 2.0; and (b) for small a (r = 0.7): a = 0.3 (full lines) and a = 0.5 (dotted lines). Straight 
lines in figures (a) and (b) represent exponential [eq. f^ ] and t^°' dependencies, respectively (in Fig. 
2a they are shown by dashed lines). 
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